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Abstract 

We present two novel dynamical features of flux tubes in the 
Friedberg-Lee model. First the fusion of two (anti-)parallel flux tubes, where 
we extract a string-string interaction potential which has a qualitative simi¬ 
larity to the nucleon-nucleon potential in the Friedberg-Lee model obtained 
by Koepf et ah Furthermore we show the dynamical breakup of flux tubes 
via gg—particle production and the disintegration into mesons. We find, as 
a shortcoming of the present realization of the model, that the full dynami¬ 
cal transport approach presented in a previous publication fails to provide the 
disintegration mechanism in the semiclassical limit. Therefore, in addition, 
we present here a molecular dynamical approach for the motion of the quarks 
and show, as a first application, the space-time development of the quarks and 
their mean-fields for Lund-type string fragmentation processes. 
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1 Introduction 


During the collision of heavy nuclei at ultra-relativistic energies it is believed that the 
energy density reached may be high enough for a quark-gluon plasma (QGP) to be 
formed. Originally such collisions were assumed to be highly transparent ii and 
the two nnclei pass throngh each other keeping abont half of their energies, thongh 
this might happen only at very high energies. On the other hand some of the indi¬ 
vidual target and projectile nncleons obtain a net color charge dne to the exchange 
of gluons in the reaction zone. As a consequence, color flnx tnbes connecting these 
color charges, are built up as the nuclei recede. If the ions have sufficient energy, 
a large nnmber of flux tubes may be produced, that may coalesce to a giant color 
flux tube or simply to an ensemble of individnal flux tubes and color ropes (highly 
charged flnx tnbes) p, ||, || . 

As kinetic energy of the ions is transferred to the color field between them it 
is believed that qnark-antiqnark and glnon pairs are created from the field via the 
Schwinger mechanism. If the pair creation rate is high enongh, the parton density 
prodnced may be sufficient for a QGP to be formed. As the QGP expands and cools, 
hadronization begins. After hadronization is complete the hadron gas expands nntil 
the strong interactions freeze ont and free hadrons are observed. 

The present status of model descriptions of this scenario is the following: On 
the one hand there are the string- and parton-cascade-models (Fritjof P, Venus 0, 
RQMD 0, HIJING p, parton cascade J^), that are based on phenomenological 
arguments and (or) pertnrbative QGD. However, all the nonperturbative mechanisms 
like string formation and decay and the final hadronization, are handled purely phe¬ 
nomenologically, i. e. by fragmentation prescriptions for the string breaking from £t- 
ting data or by coalescence models |^, |^. In order to gain more insight into 

the latter mechanism, effective models of QGD can be used, like the MIT-Bag model 
T3| , the NJL-model |]G, lAj or the various color dielectric models [IT, 17, [18 


Within these models there are attempts to describe the prodnction of qnark-antiqnark 
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pairs in flux tubes ^ or the collapse of flux tubes [Q. The main topics at¬ 
tacked in these works are the influence of the hnite transverse size of the flux tube on 
the Schwinger pair production rate, or neutralization times for the strings via pair 
production, but not the full disintegration of flux tubes into physical mesons, which 
is needed in a consistent transport description of heavy ion collisions. 

In this paper we consider flux tubes constructed in the Friedberg-Lee model. 


which has had many successes in other applications []Ig and is a fully dynamical 
held theory with conhnement. The quark degrees of freedom are handled via a 
transport equation derived by Elze and Heinz (see also Blttel et ah |2^). In 
previous publications [^, |2^ we have investigated the dynamical conhnement 

properties in nucleon-nucleon collision in the colorless version of the model as well 
as color-excitations in the full model, whereas here the outline is as follows: 

In section 2 we give a brief introduction to the Friedberg-Lee model, and show 
how the parameters of the model can be hxed by static properties of hadrons and hux 
tubes. In section 3 we extract a string-string interaction potential by a dynamical 
fusion of two (anti-)parallel strings and in section 4 we present the full transport dy¬ 
namical approach to the disintegration of a hux tube into mesons via quark antiquark 
particle production. After the discussion of some intrinsic problems connected with 
the semiclassical limit of the model, we give the space-time description of a string- 
fragmentation process of the Lund-type within a molecular dynamical approach. 


2 The Friedberg-Lee model 


The nontopological soliton model of Friedberg and Lee [^, |^, further developed 
by Goldham and Wilets has enjoyed considerable interest in recent years, 

because it provides a dynamical description of hadrons as solitons of a selhnteracting 
mean held. Static solutions of this model have been studied extensively during the 
last decade |^, |^, |^, leading to a satisfactory description of hadronic properties. 
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In its original version the Lagrangian is given by 


+ ^qa + C-a + Cq , ( 1 ) 

where the various terms are as follows: 

Cq = ^{i'-f^D^ ( 2 ) 

describes the quarks as Dirac particles of current mass mo ~ 10 MeV. The fermion 
wave functions have 4 (Dirac) times 3 (color) times n (flavor) components and the 
covariant derivative is where are the gluon helds that are 

coupled to the quarks covariantly. The term 

= ^goa'^ (3) 

describes the coupling of the quarks to the cr-£eld, which is assumed to mimic the 
long range and nonabelian effects of multi gluon exchange. The kinetic and potential 
part for the a held 

C = (4) 

contains the self interaction potential U (a) that is supposed to be quartic in the held, 

U{(^) = -^(7^ + + -^(7* + B . (5) 

The constants in (|]) are adjusted in such a way that U{a) has a minimum at a = 
0 and another, energetically lower one, at the nonzero vacuum expectation value 
(Tyac, where the potential is assumed to vanish. The constant B is chosen such that 
U{<Jyac) = 0. Since also U{0) = B, B can be identihed with the MIT-Bag constant, 
or the volume energy density of the cavity. 

Color gluon helds are introduced as in QCD except for an interaction with the 
soliton held through a dielectric function k((t), 

= -\<<7)F;,Fr , ( 6 ) 
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with fi:(0) = 1 and ^(a^ac) = 0. ^(cr) is not nniqnely defined in this model, and a 
choice must be made as to its functional form. Several suggestions have been made 
in the past, that may be summarized to the general form 


Knm{^) = \l-{^rrQ{^vac-CT). 


(7) 


Friedberg and Lee |]^ originally proposed n = m = 1, while others [^, Q have 
suggested n = 1, m = 2 or even {n,m) = (2,1), (2,3/2) and (2,2) [^. We prefer 
the parameters {n,m) = (1,2), for this choice guarantees a vanishing derivative of k 
at (Tyac' 

Color confinement is obtained by the general properties of k] it can be shown that 
a medium with vanishing vacuum value of the dielectric constant is color confining, 
i.e. net color charges would lead to infinite energy configurations ig. 


Since k and the cx-field are supposed to represent the nonperturbative structure 
of QCD, the gluon fields are treated as Abelian (one gluon exchange) Maxwell 
fields, which obey 


with 


fJ.U 


= d^Al - d^Al , 


( 8 ) 

(9) 


with the color current density 


J 


a 

V 


ya 


( 10 ) 


All the non-Abelian features of QCD are assumed to be contained in the dielectric 
function K{a). The other two field equations derived from the Lagrangian are then 

- Wv^Al) - mo - goa)'^ = 0 (11) 


for the quarks and 

B^a^a + U'{a) + \n'W)F;,Fr + goP, = 0 (12) 

for the (T—field, with the scalar density ps = TT. The primes on IF and on k' denote 
differentiation with respect to a. 
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We use the following approximations: The a and the ^“-field are treated as 
classical fields (mean field approximation). For the gluon field we choose the Coulomb 
gauge V ■ = 0 resulting in 


Vk, 


-Jo 


- - V x {kA'^ x —) = -j + nVdtA^ . 


K 


(13) 

(14) 


In the following equation (|T^ is neglected as it can be shown, that currents within a 
flux tube do not produce a magnetic field because the displacement current is exactly 
cancelled by the convection current if the string radius stays nearly constant 0. As 
a consequence, we determine the colorelectric held instantaneously by = — VAq, 
which is exact in pure 1+1-dimensional electrodynamics. 


The densities needed as sources for the equations (jT^) and (jT^) are provided by 
a transport equation derived by Elze and Heinz in the semiclassical limit of 
an exact equation of motion for the quantum Wigner function making use of the 
Dirac equation ([TI|) . The corresponding equations for the phase space distribution 


functions for quarks (/) and antiquarks (/) read 


- m*{df,m*)d^)f{x,p) = g^Pf,Ff^'^dPf{x,p) (15) 

{p^d>^ - m*{df,m*)d^)f{x,p) = -g^p^F>^‘'dPf{x,p) . (16) 

We dehne the effective mass m* = tuq + g^a, the energy oj = a/jF + and a factor 
?7 = 4 describing spin and isospin degeneracy. In terms of / and / the scalar density 
is then given by 

Ps = [ d^p—{f{x,p) + f{x,p)) , (17) 

(ztt)-^ j uj 

and the color charge density by 

P = jo = J d^p{f{x,p) - f{x,p)) . (18) 

In deriving (^) and (^) the Abelian limit of the model has been used to reduce the 
color octet components of the currents and helds to one component denoted by jo 
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and A^. Since all the nonabelian effects of QCD are assnmed to be modeled by the 
cr—field, this assnmption is jnstified. We therefore drop the color index here and in 
the following. 


The eqnations (p!5D and (pTl) are a set of nsnal Vlasov eqnations describing the 
motion of charged particles in a selfconsistently generated scalar and vector field, 
determined by eqnation (|T^ and (|T^), respectively. They are solved by the so called 
testparticle method developed by Wong et ah [Q, where every physical (anti-)qnark 
is represented by an ensemble of testparticles 


Nt 


f{x,p,t) 

2 - 

(19) 

f{x,p,t) 

Nt 

= Xi{t))6{p- Pi{t)) , 

(20) 


2=1 


with the testparticles moving according to the Hamiltonian eqnations of motion: 


dxi 

dt 

dpi 

dt 


Pi 

UJi 


m*{xi) 


iOi 


g^E{Xi) . 


( 21 ) 

( 22 ) 


The testparticle coordinates for the antiparticle ensemble obey the same eqnations of 
motion except for —g^. The simnlations are rnn in the so-called fnll ensemble 
method. 

Since we work in the abelian limit of the theory, the color charge of the qnarks 
can only be positive or negative. While this presents no problems for the mesons, 
color nentral baryons, which consist of three qnarks, cannot be constrncted in a 
straightforward way. We therefore work in the latter case with diqnarks, which carry 
baryon nnmber B = -1-2/3 and negative color charge. Thns, in the case of mesons, 
/ denotes the phase-space distribntion of antiqnarks with B = —1/3 and negative 
color charge, whereas in the case of baryons / denotes the diqnark distribntion, again 
with negative color charge, bnt B = -1-2/3. 
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2.1 Static limit 


Let us first examine the static limit of the transport equations. From quantum 
mechanics we know that the color charge density in a hadron has to vanish locally 
1 ^ 

]o{x) =< N\Q{x)\N >= J d^p{f{x,p) - f{x,p)) = 0 , (23) 

because the nucleon ground state \N > is a color singlet and Q an (arbitrary) color 
octet operator. From that we conclude that there is no colorelectric field in the 


groundstate, thus / = /. As shown by Vetter et al. |2^, the distribution functions 
have to be of a local Thomas Fermi type 


/(a:,p) = f(x,p) = , 


(24) 


where we have introduced the Fermi energy p for the quarks. After inserting the 
local Thomas Fermi distributions (|^) into the integral expression (pTl) the fermions 
can be integrated out to give 


Ps = = 

for the scalar density and 


2pm* 


TT^ 


ppf + {m*) log( 


m 


P + Pf' 


Q{p-m*) , 


(25) 


P^ 


= = 




- m') 


(26) 


for the quark density. Thus we are left with determining the soliton solution of the 
remaining cr-field equation (|I^). The equations (|^ , |^) differ from (28) in ||^ by a 
factor of 2 because of the antiparticle contribution. 

The soliton solution is constructed using a shooting method of Van Wijngaarden- 
Dekker-Brent [^. We do not intend to investigate the total parameter dependence 
of the properties of the soliton solutions, since these have been extensively discussed 
elsewhere ffl. Instead we restrict ourselves to adjust the model parameters in 


a way to reproduce the quark-number, mean mass of delta and nucleon and the 
rms radius of the nucleon for a typical baryon. With the parameters of this best 
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fit kept fixed, we use the Fermi energy fi in order to go from the baryonic solution 
{Nq = f d^rd^p(f(x, p) + f{x,p)) = 3), resembling a quark-diquark conhguration to 
the mesonic solution {Nq = 2), resembling a quark-antiquark conhguration. The 
parameter set used to obtain these solutions is displayed in table |^, where we used 
fi = 1.768for the baryon and p = 1.9582/m“^ for the meson. The physical 
properties of these soliton solutions are shown in table and where the meson 
mass, which typically is too high in these models, can be signihcantly reduced by 
the inclusion of momentum projection methods and the colormagnetic energies 0. 
The remaining properties are very well reproduced within our model approach. These 
groundstate solutions will be used in the following sections for dynamical applications. 


3 The string-string interaction potential 

Gluonic interactions in soliton bag models have been considered in two ways in the 
literature. On one hand the spin dependent one gluon exchange (OGE) interaction 
has been used to £t the strong coupling constant of the model to the experimentally 
observed N — A—splitting. In these calculations a gluonic propagator developed by 
Bickeboeller et ah [Q is used, resulting in a value for ~ 2 |0, |4T| . 

Although it is in principle possible to determine the OGE matrix elements for dy¬ 
namical systems without any special spatial symmetry, the calculation of the prop¬ 
agators is very involved and has therefore not yet been performed. On the other 
hand there have been flux tube solutions constructed with an infinite (cylindrical) 
geometry assuming a large mass for the quark-antiquark pair at the endcaps of the 
tube (Born-Oppenheimer approximation) and an instantaneous colorelectric held be¬ 
tween them 10, 01 . Within these approximations the most important property of 
the tube, the string tension r ~ 1 GeV/fm, has been htted, where again values of 
the strong coupling constant ~ 2 have been obtained. 

We proceed differently here and calculate first the response of the colorelectric 
held to an adiabatic separation of the quark-antiquark pair of the meson ground- 


9 




state, i. e. the quark and antiquark move apart with a given constant velocity which 


is taken to be Vq = —Vq = 0.1c in order to guarantee adiabacity in our calculations. 
The quark degrees of freedom are taken as frozen, so that the spatial shape of the 
quark distributions is not changed during this process. The corresponding field equa¬ 
tions 0 and (p!3D are solved selfconsistently throughout the time evolution via a 


staggered leapfrog algorithm for the complete (3-dimensional) cx-held combined 
with a (2-dimensional) hnite element method for the colorelectric held exploit¬ 
ing the cylindrical symmetry of the system. In this way we dynamically generate a 
color dielectric hux tube with a length of 8 fm and a width of ~ 1 fm as described 
0 d hx the strong coupling constant to a string constant of r = 1 GeV/fm, 


m 


obtaining = 1.92, in agreement with the earlier estimates discussed above. 

Similarly we can construct a two hux tube conhguration with either parallel or 
antiparallel colorelectric helds as schematically depicted in hgure |l|. Starting with 
this conhguration we let the quarks inside the diherent tubes move with a constant 
and hxed fusion velocity Vf = 0.1c towards each other keeping their spatial shape 
hxed througout the fusion process. We have checked that at this low fusion velocity 
the (T—held follows the quark motion without retardation. We hnally extract an 


interaction potential Vad for the two strings via 


Vad{p) = E„{p) + Egi{p) - E„{oo) - Egfoo) , (27) 

with 

Ea = J d^r + U{a) + goapsia)^ , (28) 

Egi = d^rK{a)E^ . (29) 

The interaction potential is taken to be the diherence in the held energies compared 
to inhnitely separated strings as a function of the transverse distance pt of the tubes. 

Due to the loss of the cylindrical symmetry in this scenario, we adopt the following 
prescription for determining the colorelectric held: We take the axis connecting the 
quark and antiquark of one of the hux tubes as the symmetry axis. The colorelectric 
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field is then calculated taking only the color charges within this tube into account, 
resulting in an ’undisturbed’ held conhguration Eu, as for instance shown in one of 
the tubes in the upper picture of hgure ||. The total colorelectric held is assumed to 
be given by 


E{x, y, z) = Eu{x - pt, y, z) ± Eu{x + p*, p, z) , (30) 

where the plus sign is valid for parallel tubes and and the minus sign for antiparallel 
ones. The index u denotes the ’undisturbed’ held conhguration and pt stands for the 
(transverse) distance of the symmetry axis of the individual tube to the x-axis. We 
consider this coherent addition of the undisturbed colorelectric helds to be a strong 
assumption, since we neglect the direct interaction of the individual quark with both 
quarks of the neighboring hux tube. However, in the limit pt 0, the cylindrical 
symmetry is restored and our prescription is valid again. Furthermore, the equation 
for the (j—held is solved with the full 3+1-dimensional technique developed by Vetter 
et ah m- A complete 3-dimensional hnite-element treatment of both, the Poisson 
equation (JT^) and the cr—held equation (|^), is in preparation . 


The result for the parallel conhguration is shown in hgure |^, where one can see 
that there is no long range interaction in this model due to conhnement for 
distances larger than 1 fm. But as soon as the tubes get in touch with each other 
at distances less than 1 fm, the conhnement wall between them breaks down and 
we gain surface energy from the a—held (^). The short range behaviour is strongly 
repulsive, which is clarihed in hgure where the held conhgurations are shown for 
the cr—held (left side) and the colorelectric held (right side) throughout the adiabatic 
fusion: when the charges start to overlap, the coherent addition of the colorelectric 
held leads to a doubling of the held and thus to 4 times the colorelectric energy 
in the complete overlap conhguration. This rise in energy results in the strong short 
range repulsion of the parallel strings. The adiabatic string-string potential is quali¬ 
tatively very similar to the nucleon-nucleon potential in the same model obtained by 
Koepf et al. [Q], which has been calculated as the energy diherence of a deformed 
six-quark bag and two noninteracting separated nucleons. We even get quantitative 
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agreement with those calculations if we scale the potential with the length L of the 
flux tube, which is shown in figures ^ and |^. 

For the antiparallel configuration (figure ^ again there is no long range inter¬ 
action. However, as soon as the confinement wall of the flux tubes starts to break 
down, we get an attraction coming from the decrease in energy of the a—held. An 
even larger attraction results from the colorelectric held, since the two opposite color 
charges at the end of the tubes neutralize each other, so that the color held vanishes 
at very low transverse distances and the hux tube collapses. 

Our model then supports the independent string picture, which is, for instance, 
phenomenologically used in the Lund code |^. That means that the parallel strings 
repel each other, building an ensemble of single individual strings. Antiparallel 
strings, on the other hand, do prefer to fuse and neutralize the color charges at 
the endcaps, leading to a collapse of the hux tube. 


4 Flux tube breaking 


In this section we show how a hux tube of the Friedberg-Lee model breaks up due 
to quark-antiquark pair production. One usually describes such processes with the 
Schwinger formula 0, which gives a constant pair production rate of electron- 
positron pairs in QED, depending only on the absolute value of the electric held. How¬ 
ever, this formula cannot be naively transferred to the QCD case of gg—production 
in a hux tube for very diherent reasons. 

First of all, the back reaction of the produced pairs on the external held has 
to be considered H. since this screening of the held is hnally responsible for the 
breakup of the strings. Furthermore one must consider the modihcations of the pair 
production rate due to transverse and longitudinal conhnement & EH, O, ^ ® 


Besides many other issues, these authors have shown in various bag models, that the 
production rate is strongly suppressed due to the reasons mentioned above. 

However, a detailed analysis of the time dependent pair production process within 
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the Friedberg-Lee model has not yet been performed. Therefore we have to rely on 
different phenomenological arguments, providing us with a simple guideline to the 
space-time evolution of a fragmenting string. 

The earliest and most successful model of string fragmentation is the Lund-model 
||5CI|| . Within this model one assumes the quarks to be massless and therefore mov¬ 
ing with the speed of light. Even the original QQ—pair, generating the string, is 
supposed to travel on the lightcone. Within these assumptions it is sufficient to de¬ 
scribe the fragmentation process by a 1-|-1-dimensional space-time geometry, for it is 
causally impossible that the produced gg—pairs propagate towards the endcaps (the 
QQ-paii). 

Four of the most simple string fragmentation diagrams emerging from this model 
are shown in the figures a) - |^ d). The figures show the propagation of the original 
QQ—pair and the consecutively produced gg—pairs under the influence of a constant 
and therefore confining colorelectric field. In figure ^ a) we show the breakup of the 
string due to the production of one pair in the center of the string and the symmetric 
and equal-time production of two pairs in figure b). In both cases excited meson 
states emerge, travelling with the rapidity of the orginal quarks, whereas in the sec¬ 
ond case an excited ’Yo-Yo’-state with vanishing velocity is formed in addition. In 
figure H c) and d) the breakup of a string via equal- and non-equal-time production 
of three gg—pairs is presented. In the first case the final excited ’Yo-Yo’-modes again 
have a vanishing total momentum whereas in the latter case their total momentum 
is nonvanishing, with its value depending on the difference in production time. 

There are two physical properties of the final mesons that are fixed by the model 
parameters P, the mass of the final fragments can be fixed by choosing the 
space-time breakup points appropriately, according to 

-x^ = 2mVr^ , (31) 

where m is the mass of the final meson and r the string-constant. The second property 
to be fixed is the rapidity distribution of the final fragments, that can be chosen 
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through the so called fragmentation function p, . Within the range of only a few 
parameters, that are fitted to data, a very reasonable description of particle 

spectra for hadron-hadron and heavy-ion collisions can be obtained. Motivated by 
these results we present in the following a fully dynamical space-time description of 
the fragmentation processes shown in hgure ^ within our model approach. 


4.1 The Transport Dynamical Model 


As a prerequisite we assume that the original QQ—pair travels with the speed of 
light along the ± 2 ;—axis, irrespective of the interior dynamics, which is generated 
by solving selfconsistently the transport equations (^) and (plbf ) together with the 
mean-field equations for the cr—field (0) and the colorelectric field (p!^ with the 


numerical techniques described in the previous section. The source densities enter¬ 
ing the field equations are determined via the integral relations ([T7|) and (|1^). The 
produced gg—pairs are inserted into the dynamically evolving flux-tube by assuming 
their groundstate spatial shape that has been determined in the previous section 


and with a vanishing total momentum [^. Each quark of the produced gg—pairs is 
represented by a testparticle ensemble with N-r = 50000. 

In figure § we show the color charge density p = jo, the energy density of the 
colorelectric field e = and the scalar cr—field along the longitudinal and trans¬ 

verse axis. In the upper row, at t = 1.8 fm/c, the quarks have a distance of 3.6 fm 
and the colorelectric field between them together with the confining cavity of the 
cr—field builds up. At t = 2 fm/c the gg—pair is inserted and instantaneously torn 
apart by the force of the colorelectric held. In the second row, at t = 2.6 fm/c, we 
see that the separating gg—pair causes a screening of the colorelectric held in the 
region between them, where the energy density of the color held already vanishes. 
Consequently, but with a little delay in time, the cr—held is beginning to snap oh in 
the screening region. In the lower two rows, at t = 3.4 fm/c and t = 7.0 fm/c the 
free propagation of the hnal string fragments can be seen, where we observe some 
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slight late oscillations of the a—held due to its inert behaviour caused by the large 
glueball mass. Overall, the observed time development corresponds to the expected 
space-time behaviour shown in hgure ^ a). In the lower row of hgure ^ we recognize 
a small dispersion of the q and q distributions, which is caused by the selhnteraction 
of the quarks, leading to a screening of the charges: a charge fragment at the end of 
one of the two substrings feels a smaller colorelectric force due to its interaction with 
the charge fragments in front of it. 

As a further example we show in hgure the space-time evolution of the string 
fragmentation via equal time production of three gg—pairs, in analogy to hgure 
c). In the upper row at f = 2.8 fm/c, the string is already extended to a length of 
5.6 /m; the constant colorelectric held along the string is rehected in the constant 
held energy between the charges. At t = 3.0 fm/c the three gg—pairs are inserted at 
z = 0 fm and = ±2.0 fm with a vanishing relative velocity. 

In the second row at f = 3.6 fm/c we once again see the screening ehect from 
the motion of the quarks inside the hux-tube and consecutively the beginning of the 
snapping oh of the a—held. In the following we observe, as in the previous example, 
the formation of the outer two meson pairs propagating on the lightcone. In the inner 
region of the hux-tube, the respective quark pairs penetrate each other, having al¬ 
ready acquired a small dispersion. At the time t = 4.2 fm/c, the colorelectric held is 
almost completely screened. This abrupt change of the colorelectric held, being also 
a source of the a—held (see eq. (0)), leaves the a—held highly excited. Therefore the 
cr—held snaps oh and starts to oscillate around the nonperturbative vacuum value 
cr^ac- These oscillations can be followed in the last two timesteps of hgure |^, where 
cr ~ 0 in the center of the string at z = 0 fm and a ~ a^ac around. Comparison with 
the leftmost column in hgure ^ shows that the cr—held, which is supposed to localize 
the quarks, cannot prevent a further dispersion of the color charges and the latter 
neutralize and dissolve along the z-axis to a length of about 8 fm at t = 8.4 fm/c. 
As a consequence the scalar density of the quarks is reduced by almost an order of 
magnitude. 
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In this case the a—held undergoes a complete phase transition to the nonpertur- 
bative vacuum, since this is energetically more favorable, which is clarihed in hgure 
There it can be seen, that the global minimum of the effective selhnteraction 
potential 

Ueff{,cr) = U{a)+ goaps (32) 

is at a vac, if the scalar density drops below 1/3 of the groundstate density. After re¬ 
turning to a vac, the held undergoes damped oscillations, irrespective of the remaining 
gg—fragments. 


4.2 The Molecular Dynamical Model 


As already mentioned, in the present numerical realization of our model each test- 
particle of the respective ensembles interacts not only with the testparticles of the 
other quarks, but also with the testparticles of the same quark. This selhnteraction, 
which leads to the dispersion of the quark distributions, becomes dominant in the 
limit of only a few interacting physical particles; in the extreme case, a single color 
charge distribution with an almost vanishing mass would immediately dissolve due 
to the selhnteraction. Unlike in other models, like e. g. the Walecka-model, which 
also contains an attractive scalar held and a repulsive vector held and which has 
been treated successfully in the present case the attractive a—held acts only 
in transverse direction on the produced quarks. Its longitudinal ehect is completely 
suppressed by the colorelectric held. Thus we are forced to formulate a new dynam¬ 
ical approach, that conserves the identity of the diherent quarks and antiquarks and 
is free of the selhnteraction problem. 

A possible way to overcome the problems mentioned above is to treat the quarks in 
a molecular dynamical framework. In this approach each physical particle (quark or 
antiquark) is simulated by just one testparticle and not by an ensemble of testparticles 
{Nt = Nphys) in the equations (|^ and (|20|). Due to the fact that point-like particles 
cause short range divergences in the held equations, a specihc spherically symmetric 
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distribution is assigned to each individual testparticle. In our case we choose this 
distribution to be the meson groundstate distribution, as we have calculated in the 
previous section. This provides us with the scalar density and the charge density p 
of each individual (anti-)quark. Although violating Lorentz covariance, we keep the 
spatial shape hxed throughout the time evolution, as a hrst approach to the reaction 
dynamics. A testparticle then describes the motion of the center of its distribution 
according to the Hamiltonian equation of motion (|^) . The initial momentum of 
the produced (anti-)quarks is chosen to vanish |^, and propagated according to the 


equation (^) afterwards. The total scalar density and charge density entering the 


held equations (|T^ and ([T^) is given by the sum over all (anti-)quarks. The held 
equations, on the other hand, are simulated by the same techniques as in the full 
transport dynamical model. 

In the following we show how we obtain a space-time description of the diagram 
of hgure d), thus presenting a fully dynamical description of the disintegration of 
a hux-tube into a multiple (excited) meson state within the Friedberg-Lee model. 

The corresponding space-time evolution of this process can be followed in hgures 
^ and ro, where once again the color charge density p, the energy density of the 
colorelectric held e and the a—held is presented at eight diherent time steps, which 
are chosen such, that all relevant steps of the diagram ^ d) can be followed. We start 
our presentation at f = 2.7 fm/c, shortly after the hrst gg—pair has been inserted in 
the coordinate center. We see that the quarks of this hrst gg—pair already propagate 
along their respective directions and observe a prompt screening of the colorelectric 
held in the center, whereas the more inert a—held has not yet changed signihcantly. 
In the second row of hgure ^ at f = 3.2 fm/c, the remaining two gg—pairs have 
just been inserted at = ±2/m at f = 3.0 fm/c and are moving under the inhu- 
ence of the colorelectric held. We see how the screening proceeds not only at the 
coordinate center, but also between the other fragments. The conhning a—held still 
shows no reaction to the interior dynamics. In the third row we clearly see how the 
two lightcone fragments have been formed as well as the two excited ’Yo-Yo’-modes 
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in between after bAfm/c. The motion of the latter at this particnlar timestep is 
snch that the distribntions almost completely sit on top of each other, so that the 
colorelectric held vanishes. 

We additionally clearly see how the conhning a—held starts to form aronnd the 
hnal fragments. This process is completed after 7.2 fm/c as we hnd in the lower row 
of hgnre P, where the a—held again undergoes slightly damped oscillations in the 
regions of the nonperturbative vacuum. Due to the diherent production times of the 
inner gg—pairs, the respective pairs of the hnal fragments have a nonvanishing total 
momentum when forming an excited meson and thus start to propagate along the 
z-axis as can also be seen in this hgnre. 

The late time evolution of this reaction is shown in the hgnre 0 at the time steps 
t = 7.8, 8.6,10.6,12.2/m/c; note the rescaled z-axis for presentational purposes. 
During this late time development we hnd a more or less complete correspondence 
of the dynamical behaviour of our model with the 1-1-1-dimensional Lund-diagram 
^ d). We see the propagation of the hnal fragments as well as the remaining ’Yo- 
Yo’-modes of the inner two excited mesons. The latter are obviously undamped as 
expected from the simple diagramatic representation. 


5 Discussion, Summary and Outlook 


We have presented two novel features of hux tubes in the Friedberg-Lee model. First 
we have extracted a string-string interaction potential by an adiabatic fusion of two 
parallel or antiparallel strings. The potential is obtained by comparing the colorelec¬ 
tric and the a—held energies of the fused conhguration to the (inhnitely) separated 
conhguration. For the parallel strings it shows a qualitative similarity to the nucleon- 


nucleon-potential of Koepf et ah which has been calculated in a similar way as 
the energy diherence of a deformed six-quark bag and two noninteracting seperated 
nucleons. Therefore it seems to be a general feature of the Friedberg-Lee model, 
that the interaction of two color neutral objects (nucleons, parallel strings) can be 
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described in the following way: 

Due to confinement in this model there is no long range interaction (over dis¬ 
tances typically of the order of 1 fm and more). In the intermediate range, when the 
confinement wall starts to break down, we have a balance of gaining surface energy 
from the a—field and increasing the total energy from the coherent addition of the 
colorelectric field. In the short range region the potential is strongly repulsive due 
to the unscreened quark color charges at the endcaps in the case of parallel strings, 
and in the NN-interaction due to the unscreened (aligned) color spins of the quarks. 
In the case of antiparallel strings, however, a strong attraction is found leading to a 
neutralization of the color charges and the collapse of the flux tube. 

Secondly we have shown how the flux-tube breaking proceeds in the Friedberg- 
Lee model. Since no estimates of the pair production rate in this model have been 
performed so far, we give a space-time description of some of the simplest Lund-type 
diagrams. We have shown the full transport dynamical behaviour of the quark and 
antiquark distributions for a string fragmenting via the production of one gg—pair 
and three gg—pairs at an equal time. In the first case the space-time behaviour of 
two meson pairs propagating on the light-cone can be reproduced by our model. In 
the latter case of three produced gg—pairs, the outer meson pairs are also formed, 
but in the center of the fragmenting string, where the two parallel ’Yo-Yo’-modes are 
supposed to build up, we observe a different behaviour: 

When the respective quark distributions start to penetrate and neutralize, a small 
dispersion of these is already observed. At this time the colorelectric field is com¬ 
pletely screened in the center of the string. This sudden change of the source of 
the flux-tube leaves the cr—field highly excited and thus it starts to oscillate around 
the nonperturbative vacuum value. The a—field is not able to prevent the further 
dispersion of the quark distributions, which neutralize and dissolve along the string- 
axis. With the scalar density being significantly reduced in this case, the a—field 
undergoes a complete phase transition to the nonperturbative vacuum. 

Therefore we have to face the following problems within this realization of the 
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dynamical model: 

First of all the large value of the strong coupling constant a* = 1.92, being re¬ 
sponsible for the large energy density of the colorelectric held, leads to a very deep 
hux-tube, with values of the cr—held up to —cr^ac, whereas the perturbative region 
is thought to be around a ~ 0. Therefore the a—held is highly excited when the 
quark charges are screened, causing large amplitude oscillations that are only slowly 
damped. In order to avoid these oscillations and provide a faster damping, one might 
think of an additional coupling of the a—held to its chiral partner, the pion-held, in 


the chirally symmetric 0(4)-version of the model [^, |^, ^|. By this extension, the 
oscillations of the a—held can be used to excite low mass pion-modes, that disperse 
the energy of the cr—held ehectively. 

Additionally one could determine a new parameter set of the model. A lowering 
of the bag-constant, for instance, would facilitate the string-formation, as well as the 
fragmentation. However, the present value of B = bQMeV/fm? is already in the 
lower region of the ones that are generally used. A larger value of the scalar coupling 
constant go would facilitate the formation of bags and thus the localization of the 
quarks in a perturbative region. 

The second problem is the dispersion of the quark distributions, which is caused 
by their selhnteraction. This dispersion is an artefact of our numerical method, for 
we deal with only a few quarks that are treated as classical charge distributions and 
interacting strongly via a classical conhning potential. In this case the selhnterac¬ 
tion of the color charges becomes dominant and generates a screening of the charges, 
which dually is responsible for the dispersion of these. 

We have shown a possible way to overcome these problems by treating the quarks 
and antiquarks as point particles with an appropriately chosen charge distribution as 
commonly used in molecular dynamical simulations. We have shown within a hrst 
and simple approach, where we keep the spatial shapes of these distributions hxed, 
how the Lund-type dynamics of the string-fragmentation can be reproduced for even 
more complex Lund-diagrams. An improvement of the molecular dynamical approach 
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can be achieved by implementing a time-dependent radius or dipole-moment of the 
quark distribution. These multipole moments could account for Lorentz contraction, 
polarization effects of the distributions or monopole oscillations. 

Summarizing we can say that our model provides us with a useful tool to de¬ 
scribe the full dynamical evolution of string formation and decay via multiple quark- 
antiquark production. The aim for future investigations should be to apply our model 
to more realistic scenarios, like multiple string production and decay in relativistic 
heavy-ion collisions, for describing the formation of hadrons out of an expanding and 
cooling quark-gluon plasma. 
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Figure 1: Parallel and antiparallel flux tube conflgurations drawn schematically. 


Figure 2: The cavity of the a—held and the colorelectric held throughout the adi¬ 
abatic fusion process of parallel flux-tubes. The transverse distance of the tubes 
is (from top to bottom): p = 1.75,1.25,0.75,0.25/m. The equidistant contour 
lines start at —0.125to 0.25 fm~^ (in the lower row from —0.175to 
0.25 

Figure 3: Adiabatic string-string potential for the parallel configuration. The solid 
line shows the sum of the contributions from the a— and the colorelectric held. 


Figure 4: Adiabatic string-string potential for the antiparallel configuration. The 
solid line shows the sum of the contributions from the a— and the colorelectric held. 


Figure 5: Space-time representation of some of the simplest Lund stringfragmentation 
graphs. 


Figure 6: The figure shows the time evolution of the color charge density p, the 
energy density of the colorelectric held e and the a—held for the breakup of the 
string via the production of one gg—pair. For calculating the motion of the quarks 
the full transport dynamics have been used. The equidistant contour lines run from 
—0.5 fm~^ to 0.5 fm~^ for p, from 0.0 fm~^ to 4.0 fm~^ for e and from —0.15 fm~^ 
to 0.25 fm~^ for a. The respective charges of the generating QQ—pair are indicated 
by a -|- or — sign. 
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Figure 7: The figure shows the time evolution of the color charge density, the energy 
density of the colorelectric field and the a—field for the breakup of the string via the 
production of three gg—pairs. For calculating the motion of the quarks the full trans¬ 
port dynamics has been used. The equidistant contour lines run from —0.5 fm~^ to 
0.5 fm~^ for p, from 0.0 to 3.0 for e and from —0.125 fm~^ to 0.25 fm~^ 
for a. The respective charges of the generating QQ—pair are indicated by a -|- or — 
sign. 


Figure 8: The effective selfinteraction potential for different values of the scalar 
density, starting from the solid line to the dotted line: 0.110/m“^, 0.072/m“^, 
0.050/m-3, 0.043/m-3, 0.036/m-^, 0.021/m'^, 0.000/m-T 


Figure 9: The figure shows the early time evolution of the color charge density, the 
energy density of the colorelectric field and the a—field for the breakup of the string 
via the production of three gg—pairs at different time steps t = 2.7, 3.2, 5.4, 7.2 fm/c. 
For calculating the motion of the quarks the molecular dynamical approach has been 
used. The equidistant contour lines run from —0.35 fm~^ to 0.35 fm~^ for p, from 
0.0 fm~‘^ to 3.0 fm~'^ for e and from —0.15 fm~^ to 0.25 for a. 


Figure 10: The figure shows the late time evolution of the color charge density, the en¬ 
ergy density of the colorelectric field and the a—field for the breakup of the string via 
the production of three gg—pairs at different time steps t = 7.8, 8.6,10.6,12.2 fm/c. 
For calculating the motion of the quarks the molecular dynamical approach has been 
used. The equidistant contour lines run from —0.35 fm~^ to 0.35 fm~^ for p, from 
0.0 fm~^ to 3.0 for e and from —0.15 fm~^ to 0.25 for a. 
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parameter set 


a [/m-2] 

0.0 

b [fm~^] 

-419.3 

c[l] 

4973.0 

B 

0.283 

9o [1] 

8.0 

mo [fm~^] 

0.025 


1.768 / 1.9582 

7 ] 

4 

(Tvac 

0.253 

glueball mass [Geh^] 

1.045 


Table 1: The table shows the parameters used to construct a mesonic and a baryonic 
soliton solution of the sigma held equation. The hrst value of /i is used for the 
baryonic solution, the second one for the mesonic. 
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Parameter set 

meson 

Experimental Data 

E [MeV] 

798 

465 

RMS [fm] 

0.653 

0.66 


Table 2: The table shows the results of the hts for the meson compared to the 
experimental data. The experimental data are the mean values of the pion and the 
rho-meson (taken from |]I6| ). 


Parameter set 

baryon 

Experimental Data 

E [MeV] 

1099 

1087 

RMS [fm] 

0.693 

0.83 


Table 3: The table shows the results of the £t for the baryon compared to the 


experimental data. The experimental data are taken from [|^ 
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